On exact solutions to the scalar field equations in standard cosmology 
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A simple algebraic method to obtain exact solutions to the scalar field equations in spatially flat 
FRW cosmology is derived. The field potential fuction is reduced to two terms which can be used to 
determine some characteristic inflationary variables. To demonstrate the method, exact solutions 
resembling typical potentials of general interest are worked out. Finally, it is shown that this method 
of solution is also applicable to a perfect fluid in the same metric, for a curvature k > and fixed 7. 
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I. INTRODUCTION 

It is a common issue in Cosmology to make use of 
scalar fields <f> as the responsible agents of some of the 
most intriguing aspects of our universe (see [1] and ref- 
erences therein). Just to mention a few, we find that 
scalar fields are used as the inflaton which seeds the pri- 
mordial perturbations for structure formation during an 
early inflationary epoch; as the cold dark matter can- 
didate responsible for the formation of the actual cos- 
mological structure and as the dark energy component 
which seems to be driving the current accelerated expan- 
sion of the universe. The key feature for such flexibility 
of scalar fields (spin-0 bosons) is the freedom one has to 
propose a scalar potential V (</>), which encodes in itself 
the (non gravitational) self-interactions among the scalar 
particles. Recently, scalar fields coupled to gravity (in a 
FRW background) have also appeared in connection to 
the so called string theory landscape [2] , where the scalar 
potential V((j>) is usually thought of as having many val- 
leys, which represent the different vacua solutions. The 
hope is that the statistics of these vacua could explain, 
for example, the smallness of the cosmological constant 
(the simplest candidate for dark energy). 

This coupling of gravity with scalar fields is what is 
known as scalar field cosmology. In order to construct 
this model we start with the line element for a ho- 
mogeneous, flat and isotropic universe, the Friedmann- 
Robertson- Walker (FRW) metric 



ds 2 = -dt 2 + e 2Q(t) [dr 2 + r 2 dn 2 ] 



(1) 



where a(t) = e a w is the scale factor. The classical evo- 
lution is obtained by solving Einstein's field equations. 
Using the energy momentum tensor of a scalar field we 
get the Einsten-Klein-Gordon equations 



3H4> = - 



dV{<j>) 



(2) 



'Electronic address: marco@fisica.ugto.mx 



where H = a/a. Also, H satisfies that H=—(f> 2 /2. As 
stated at the beginning, by choosing appropriate poten- 
tials, different aspects of the universe can be modeled. In 
the following we shall use 87rG=l. 

Because of the non linear character of Einstein's equa- 
tions, exact solutions to the scalar cosmology equations 
can only be found using simple methods for very partic- 
ular potentials (constant and exponential) . But for more 
appealing potentials, that explain some of the properties 
of our universe, exact solutions can not be found. In 
order to find exact solutions for complicated potentials 
one must use elaborate techniques, like Chimento and 
Jakubi's solution for a perfect fluid [3]. However, the 
need for an exact solution to the problem has been sur- 
passed by the fact that what matters most in the model is 
the dynamical aspects of the field, and therefore, the ap- 
proximations made using the slow roll method have pre- 
vailed over the need to find such exact solutions. This ap- 
proximation essentially means that the dynamical terms 
of the scalar field, denoted by ^<fi 2 and by <f) in eqs.(2), 
are neglected with respect to the other terms in the equa- 
tions. This permits us to write down the physical quanti- 
ties in terms of the slow roll parameters e and 77. Also the 
use of numerical in connection with dynamical systems 
is common to complete the solution [4]. 

Despite the success of the slow roll approximation, it 
is clear that having a more general method to find exact 
solutions to eqs.(2) may allow us to make a better de- 
scription of the phenomenology of this particular model, 
and would give a direct way to probe the exactness of the 
numerical algorithms extensively used. 

The present work is in this direction, to develop a sim- 
ple algebraic method to find exact solutions for a wide 
variety of scalar field potentials. The paper is arranged 
as follows: in Section II the method is presented for an 
arbitrary scalar field potential, followed by a disscusion 
on dynamics in Section III. A few examples of general 
interest are presented in Section IV, and the connection 
to the exact solutions of scalar field equations in perfect 
fluid is made in Section V. Conclusions and outlook are 
given in Section VI. 
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II. THE ALGEBRAIC METHOD 



III. DYNAMICAL VARIABLES 



The system of coupled eqs.(2) has three unknowns, 
4>(t), V(4>) and a(t). To solve it, one of these variables 
has to be given a priori. It is customary to look for the 
solution for a given V{<p), but as it is known, very few 
potentials arc known to lead to exact solutions. 

In this work, it is proposed to reduce the equations to 
a simpler form first, which helps to solve them exactly. 
In order to do this, let us consider the function V a ((f>) 
defined as 



V a {<l>) = V{<j>) + \ft 



(3) 



Now, with the change of variable dt — d(f>/<p 7 one obtains 

(4) 



dV a _ dV 
d(f> d(f> 



Hence, eqs.(2) can now be rewritten as 

2 

= V a 



dVa 

' d(j) 



(5) 
(6) 



These equations are now easily solved, and the function 
V a {4>) (> 0) shows up as the main part of the potential 
function, driving the dynamics of the scale factor. To 
solve them, note that eq.(5) defines a/ a as a function of 
</>, H '((/>), which when inserted into eq.(6), gives the scalar 
field 4>{t) as a function of t, at least in quadratures 



(7) 



Finally, inserting <fi(t) into eqs.(3) and (5) gives V((j>) 
and a(t), respectively, and the solution is completed. 

Obviously, one could simply have begun by giving 
H=H((f>), but it is usually desired to have some descrip- 
tion of the potential instead, and for this reason it is 
preferable to give V a (4>). One could also use H(t) to de- 
termine <j)(t), since 



^ 2 
2 V 



dH(t) 
~dT 



implies that 



A<j>(t) = ± J ^-2 



dH(t) 



dt 



dt 



(8) 



(9) 



and since V a = 3H 2 (t ) , a complete knowledge of H(t) fully 
determines the solution to the problem. For example, if 
H(t) = a/t, then V{4>) = (3a 2 - a) exp[-2A0/V2a] and 
Acj)(t) = \/2a\nt. The difference when we proceed this 
way is that we lose control of the form of the potential 
V((f), and the solutions, for example, may not comply 
with slow roll. [5] 



The potential function V((j)) in eq.(3) consists now of 
two positively defined competing terms, V a and </> 2 /2 (see 
the solution in ref.[6]), the second term being neglected 
in the slow roll approximation. In the limit <ft — > 0, wc 
find that 



lim V{4>) = V a (4>) 

0^0 



(10) 



and the scalar field potential is just V a . 

Finding two competing terms to drive the dynamics 
of the system has a quantum mechanics equivalent. For 
the phase space dynamics of the Schrodinger equation in 
one dimension, in terms of the wave function ip{x) and 
its derivative with respect to x, ip'(x) [7], the potential 
function V(x) appears as the driving term of the dynam- 
ics, competing with the decelerating term defined by the 
energy E of the system, and leading to the bifurcation 
trajectories in phase space [7]. Moreover, the similarities 
of the scalar field equations with a non-linear Schrodinger 
problem can be further realized [8] . 

One can use V a to determine some characteristics of 
the dynamics of the system. For example, from eq.(5) 
we have that 



a 
a 



± 



and from H—d/a — (a/a) 2 =— <fi 2 /2, that 



a 
a 



l Va ~l l 



Hence d has the same sign as o when 0^0. 

To end this section, let us determine the parameters 
e and rj used in the literature to quantify the inflation 
regime. Here, wc can define them as functions of the 
potential terms 



M - d 1 



dt H 



V a 



v(t) = 



1 H 
2HH 



V a 



V 



(11) 



(12) 



2 , and the primes 
p. Hence, e and rj 



where, to simplify, wc wrote for ^ 
denote differentiation with respect to 
are the relative strenghts of the two potential terms, and 
of their derivatives. In the literature, a value of e = 2 
defines the radiation era, while e < 1 defines an infla- 
tionary epoch [9]. Therefore, V a is the dominant term of 
the potential function during inflation, a condition that 
clarifies its appearance in the solution. 



IV. FEW EXAMPLES OF GENERAL INTEREST 

Now, let us consider some examples of general interest. 
Keeping in mind that V a is the main part of the solution, 
the examples are labeled accordingly. 
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A. V a = A0 2 



D. 



V a = V e- a * 



The X(j) 2 potential is one of the basic functions used in 
the literature [10], and it is a basic potential in standard 
cosmology [11]. 

In this case, 7J 2 =A0 2 /3, and dV a /d(f>=2\(f>. Therefore, 
from eq.(7) it is easy to find that 



(13) 

0) = 0, and 
(14) 



Hence, <f> is constant. Letting (j) = (j>(to 
using eqs.(3) and (5) we get 

V(4>) = \<t> 2 - 2^ 



One of the few problems exactly solved in the litera- 
ture, is that of an exponential potential function [15]. In 
this case, the solution found with the method described 
above gives, for </> = 0, 



<j>(t) = - In 
a 



v{4>) = v \i 




a 



-OL<f> 



(22) 



(23) 



a(t) = aoe 



(15) 



Obviously, one would be tempted to pick A < in order 
to make a(t) a growing function of i, but that would make 
4>(t) an imaginary function of t. 



B. V a = A</> 4 



/ a 2 V y 
a(t)=a 1 + yY 



(24) 



The reader can recognize now that this is the solution 
used in the literature, where the two terms of V{<j>) are 
grouped together, leaving the essence of the competing 
term V a out of the picture. 



The A</> 4 potential is also one of the basic functions 
used in the literature [12], and in cosmological models 
[13]. The solution is 



<f>(t) = o e ±4 ^ 



V{4>) = \4> 4 - -Xcf, 2 



(16) 
(17) 



a(tHa exp j-|> e ±8 \/i(*-Mj ( 18 ) 

Notice the appearance of a double exponential in a(t). 
In this case, the radiation era happens when <j> = 2, and 
inflation comes when <p>2\/2. 



E. V a = Va (cOSh(/30) - 1) 

Another potential used in the literature is that of a 
hyperbolic cosine [16]. The solution found here, with 
</>o=0, is 



(j){t) = -^arcsinh (tan ( Yy/^j j ( 25 ) 



V{4>) = V (cosh(/30) - 1) - V ^- (cosh(/30) + 1) (26) 



C. Va = \<t> 2n , n > 2 



a(t) = a cos^ f y y/3 2 t 



(27) 



The cases with integer n > 2, all have the same kind of 
solution. Although this is not a very common potential 
in the literature [14], it is solved here to compare with 
the two previous cases. 



In this case, e(t) = (f3 2 /2) coth 2 (^/2). 



F. A Morse type potential 



<j>(t) 



b 2 - n ±2n(n-2)^ (t - to) 



9n 2 



a(t) = ao exp 



An 



j2-n 
b 



±2n{n-2)\j-{t-t ) 



(19) 
(20) 



(21) 



The previous examples were related to potentials used 
in the literature. They may, however, not be able to com- 
ply with the desired characteristics of the model. Nev- 
ertheless, one can now generate a variety of solutions to 
eqs.(2) that may be more useful for that purpose. 

If, for example, one would like to find a potential func- 
tion as smooth as the Morse potential V(x) = A(e 2ax — 
2e ax ) with the method described here, one can begin by 
letting 



V a 



a 2 (e^ 



1 



4 



Then, we will have that 



and 



1 _ ) e 2» _ 2e ^ + 1 



V(</>) = a 2 



By requiring, for example, that [3 2 = |, we end up with 



the potential 

V{4>) = a 2 



i e Vs<t> _ 2e ^<t>/ 2 + 1 



(28) 



This potential has its minimum at <fi m = In 2 with 

a value of V((/> m ) = —a 2 , and with V"(4> m ) > 0. It is 
sketched in Fig.(l). The value for </>^— oo is a 2 . Thus, 
the smaller the a the smoother the curve. 




FIG. 1: A Morse type potential, the kind that may satisfy 
the slow roll conditions. 

This kind of potential may be more capable of satisfy- 
ing the slow roll conditions. Moreover, by requiring the 
lower sign in the previous equations, we find that a(t) is 

a(t)=a e^(l- 2 -^) (29) 
therefore allowing inflation. 



SOLUTIONS FOR A PERFECT FLUID 
WITH k > AND FIXED 7 



not the case, that is very difficult to achieve. However, 
one can still find some cases where the method of Section 
II may be useful. 

Let us begin by noting that since 



^[iH 2 - l -4> 2 - V{<f>) + const)) = 
then any pair of equations of the form 



(30) 
(31) 



3Hcf) 



dV{(j>) 







3H 2 



\^) 2 + V(<t>) + F(t) + const. 



will have the same solutions as the spatially flat case, 
provided that F(t) = 0. We shall denote the integration 
constant appearing in these equations as A [18]. 

If we had included A in eqs.(2), we would have found 
that in order to have a potential V((f>)=\(f> 2 , A=— 2A/3, 
and that a value of A = 16A/9 would have been found 
for the potential V{4>) = \(4> 2 - 4/3) 2 . 

Now, the scalar field equations for a perfect fluid in the 
FRW metric, as used by Chimento and Jakubi [3], have 
the structure given above. Those equations are 



3tf 2 







k 



(32) 
(33) 



where the energy density and pressure of the fluid are 
given by p = po/a 3j and p = (7 — l)p, respectively, with 
< 7 < 2. 

Hence, we have that any solution to eqs.(2) is a solution 
for the pair of coupled eqs. (32-33), provided that 



k > , 
p{t) - 3 



k 



7 = 



Notice that these conditions also satisfy that 

1 12 1 k 1 12 

H = — d> 2 - -7P+ ^7 = --d> 2 
2 V 2 " a 2 2 V 



(34) 
(35) 
(36) 

(37) 



and the connection between solutions is complete. 

In the literature, a value of 7 = 2/3 in cosmology is 
related to a static universe. Here, this value appears even 
in solutions where H(t) is not constant. 



The case considered here of a scalar field in a spatially 
flat FRW metric is the simplest one in the model. If the 
equations were linear, one would expect to infer from this 
method solutions to more complex systems. Since this is 



VI. CONCLUSIONS AND OUTLOOK 

In this work, the solution to eqs. (2) proceeded by giv- 
ing an ansatz for part of the potential function, which 
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was then related to the parameters that define the infla- 
tion regime. One can also start by proposing the field 
4>{t) and from there obtain the whole solution. 

If we proposed <f> as a polynomial function of t, cf) — 
a nt n ,[l7] then we would find that the highest exponent 
would give 4> oc f™-\ then H oc i 2 ™" 1 , and V(t) oc t 4n ~ 2 . 
If we worked out the case n = 1, we would arrive to the 
solution given in Section IV. A, with the leading term of 
V{4>) proportional to <p 2 . However, for n >• 1, we would 
find that V <~ </> 4 , therefore leading to the roots of the 
\(f> 2 and A</> 4 models. 

In conclusion, we have given here an algebraic method 



to find a large variety of solutions to the scalar field equa- 
tions in spatially flat FRW metric, which may be useful in 
more complex systems. Hopefully, it will give a deeper 
insight into the phenomenological aspects of the scalar 
field model. 
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